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Abstract 

We calculate the reduced density matrix for the inflaton field in a model of chaotic 
inflation by tracing out degrees of freedom corresponding to various bosonic fields. 
We find a qualitatively new contribution to the density matrix given by the Euclidean 
effective action of quantum fields. We regularise the ultraviolet divergences in the 
decoherence factor. Dimensional regularisation is shown to violate the consistency 
conditions for a density matrix as a bounded operator. A physically motivated con- 
formal redefinition of the environmental fields leads to well-defined expressions. They 
show that due to bosonic fields the Universe acquires classical properties near the onset 
of inflation. 
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1. Introduction 



The combination of general relativity with particle physics can yield viable models for 
the early Universe. This leads, in particular, to the idea that the Universe underwent a 
period of accelerated expansion ("inflation") at an early stage. Inflation not only looses 
the dependence on initial conditions, it also provides a quantitative scenario for structure 
formation. The interesting question is of course why and how inflation arises in the first 
place. It is generally assumed that a satisfactory answer can only be found in the realm of 
quantum cosmology. In this respect, ideas about the quantum nucleation of the Universe 
from "nothing" could yield a quantitative scenario for the emergence of classical inflation. 
The problem of classical initial conditions is then replaced by the question for an appropriate 
choice of the quantum cosmological wave function. Well-known choices are the no-boundary 
proposal |TJ and the tunneling proposal 0, which are in fact a whole class of proposals. There 
exist other interesting proposals such as the SIC put forward in || , that is in particular well 
equipped to deal with the case of a classically recollapsing universe. 

Inflation is assumed to take place near the GUT mass scale that is thought to lie about 
five orders of magnitude below the Planck scale. For this reason it is not necessary to rely on 
the most fundamental level of quantum gravity (such as superstring theory), and effective 
theories such as canonical quantum gravity should give an excellent approximation. Using 
either the Euclidean path integral or the Wheeler-DeWitt equation, the "transition" from the 
classically forbidden region to the classical inflationary regime has been investigated in great 
detail. Using normalised wave functions, a probability peak for the mass scale of inflation 
has been calculated in various models 0, ||. Such a probability peak has been interpreted 
- assuming that we inhabit a generic Universe - as providing a criterium to select amongst 
the members of an "ensemble" of classical universes. 

However, quantum theory does not yield a classical ensemble. Since interference effects 
can play a crucial role, it would be inconsistent to interpret these results directly as giving 
probabilities for inflationary universes with different Hubble parameters. A complete analysis 
would thus have to include a quantitative discussion of the quantum-to-classical transition. 

It is now generally accepted that the classical properties for a subsystem arise from the 
irreversible interaction of this quantum system with its natural environment. Starting with 
the pioneering work of Zeh in the seventies 0, this concept of decoherence has been developed 
extensively, see [0] and [§] for reviews. Quite recently, this continuous loss of coherence was 
observed in quantum-optical experiments |5|. 

Quantum-to-classical transition through decoherence has also proven fruitful in quantum 
cosmology. The first application of this concept to gravitational systems was performed in 
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a Newtonian framework where the (formally quantised) gravitational field acquired classical 
properties through interaction with masses In quantum cosmology, it was suggested 

to consider background degrees of freedom (such as the scale factor) as the relevant system 
and higher-order perturbations (such as gravitational waves or density fluctuations) as an 
environment [TT|. A quantitative discussion of this idea was first done in |T!|, where it was 
demonstrated how the scale factor and the inflaton field can acquire classical properties. This 
idea was further pursued in many papers @. The quantum-to-classical transition through 
decoherence plays also a crucial role for the primordial fluctuations that eventually serve as 
classical seeds for galaxies and clusters of galaxies |L3] . 



One big problem that remained in most of the above papers was the issue of regularisation. 
Since there are infinitely many environmental modes involved, there arise formal divergences 
that have to be dealt with. This was done in most cases only heuristically through the 



choice of a more or less appropriate cutoff in the number of modes. It was argued in [14 



that dimensional regularisation is only applicable for the phase part of the decoherence 
factor, not its absolute value. This observation will be confirmed in detail in our paper. 
Motivated by the finiteness of the decoherence factor in QED, a redefinition of the fields 
was performed in a particular quantum cosmo logical model in to find a finite result. (A 
similar proposal was discussed recently in Hl6| .) In the present paper we shall find in a wide 
class of models that there exists a distinguished redefinition that renders the decoherence 
factor finite. 

A redefinition of environmental degrees of freedom changes, however, the reduced density 
matrix for the system (in particular, the system momentum is changed). This was discussed 
in ||17|| , where it was suggested to use a redefinition that eliminates that part of the decoher- 
ence effect that solely arises from the change in three-volume, so that only decoherence due 
to "particle creation" remains. It is thus clear that a physical principle has to be invoked to 
select the appropriate choice of environmental variables. 

The purpose of our paper is to explore the possibilities to get a physically reasonable, 
finite, decoherence factor that can sensibly be calculated to understand the onset of infla- 
tion. For this purpose, the non-diagonal elements in the reduced density matrix have to be 
discussed. Using the criterium of normalisability for the wave function in one-loop quantum 



cosmology, the diagonal part was already discussed in the earlier papers [18, 0, [L9L pL |20| 



where the focus was on the derivation of probabilities. The same criterium will be employed 
here for the study of decoherence. We shall apply our discussion to a variety of models. The 



case of fermions, however, is discussed in a separate paper [^T], since a new formalism has 
to be applied and novel features appear. 

Our paper is organised as follows. In section 2 we briefly review the form of the quantum 
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cosmological wave function when the higher-order modes are in their vacuum state. We 
then give a derivation of the correct normalisation for this state, which is different from 
a (wrong) normalisation that is sometimes used in the literature. In section 3 we discuss 
the diagonal part of the reduced density matrix for the inflaton field and, in particular, 
its dependence on the one-loop effective action. Sections 4-6 comprise the main part of 
our paper. In section 4 we discuss general properties of the non-diagonal elements of the 
reduced density matrix. In section 5 we demonstrate that the application of dimensional 
regularisation violates basic properties of the density matrix. Section 6 then shows that there 
exists a distinguished redefinition of the environmental degrees of freedom that renders the 
decoherence factor finite. Section 7 is devoted to a physical interpretation of the results 
obtained and an outlook on future work. 

2. Cosmological wave function 

We consider the quantum-cosmological wave function in the context of chaotic inflation. 
This includes all cases of boundary conditions mentioned in the Introduction [l], [|, [3| . The 
background variables are the scale factor a of the inflationary universe and the inflaton field 
(p. In addition we assume the presence of spatially inhomogeneous fields /(x) that may be 
the higher modes of the inflaton or some other field. This fields play the role of microscopic 
(environmental) modes and are treated perturbatively. They will thus be integrated out 
below to study their decohering influence on the background variables. 

We shall consider the cosmological quantum state in the framework of the reduced phase 
space quantisation, that is, in the representation of physical variables. This has the advantage 
that a well-defined inner product is available for the interpretation of the density matrix. We 
must emphasise, though, that this reduced quantisation is limited to semiclassical branches 
of the wave function, but this is the framework where the present investigation takes place. 

In this representation the role of /(x) is played by physical polarisations of linear fields 
of all possible spins, and there is one physical degree of freedom in the homogeneous sector 
of (a,(p), which without loss of generality can be identified with (p. Thus, the full set of 
physical variables is 



With the decomposition of f(t, x) into a discrete series of spatial ortho normal harmonics 
Q n (x) on a section of a closed three-space, 




(2.1) 




(2.2) 
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<j>(x) can equally be represented by the countable set of (ip(t), f n (t)). 

Other variables (including lapse and shift functions, unphysical components of linear 
vector and tensor fields) are parametrised in terms of physical ones within a particular 
choice of gauge conditions fixing the local gauge and coordinate symmetries. In the sector 
of (a, ip), for example, in the cosmic-time gauge this parametrisation can have a very simple 
form in the slow-roll approximation, corresponding to slowly varying (practically constant) 
inflaton field. It can be taken to coincide with an approximate classical solution with the 
cosmic time t, 

<p(t) = <p, (2.3) 

a(t) = — cosh Ht, H = H(cp), (2.4) 
H 

where H(ip) = 87rV(<p)/3m 2 P is the Hubble constant generated by the inflaton potential V((p). 
We emphasise that the time parameter that appears in the reduced formalism is equivalent 
to the standard WKB-time parameters that appears in the semiclassical approximation to 



the Wheeler-DeWitt equation |22 . 



In the no-boundary or tunneling prescription of the cosmological wave function, the 
classical background fl2.4| ) arises as an analytic continuation of the solution of Euclidean 
Einstein equations, 

<p(r) = <p, 

a(r) = — sin Ht, H = H(<p), (2.5) 
H 

into the complex plane of Euclidean time r, 

r=^ + U. (2.6) 

This is usually interpreted as a "quantum nucleation" of the Lorentzian DeSitter spacetime 
from the gravitational instanton - the Euclidean four-dimensional hemisphere with radius 
R = 1/H((p). In more standard language, the Euclidean section of this instanton just cor- 
responds to a classically forbidden region, and the analytic continuation into the Lorentzian 
regime corresponds to the emergence of time from timeless quantum gravity p3| . The cosmo- 
logical wave function in Lorentzian time can also be understood as an analytic continuation 
of the Euclidean wave function, which in the semiclassical one-loop approximation has the 
form 

Mr\<P,f) = ^=e- I ^Y[^(r^\f n ) . (2.7) 
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Here I(r,ip) denotes the Euclidean Hamilton- Jacobi function for the solution (|2.5| ), and 

€{r, <p\fn) = -4= exp (-la%r)^\fA , (2.8) 

where u n (r) is the set of basis functions of the Euclidean linearised equations of motion for 
the Euclidean action I[<p, f] of physical variables; it satisfies the regularity condition on the 
Euclidean hemisphere containing the pole r = 0, 

F n (d/dr)u n (r) = 0, u n {r) = reg, (2.9) 

F n (d/dr)6(r - t') = .f/^/j ■ (2.10) 

Analogously, u v (t) in (|2.7|) denotes the basis function of a spatially homogeneous inflaton 
mode. This expression was obtained both by solving the Wheeler-DeWitt equation and by 
calculating the Euclidean path integral |19|, |4], |25], |26|]P]- It is important that in the above 
expressions the basis functions and Gaussian states of quantum variables f n are all calculated 
on the background of the classical solution parametrised by the inflaton variable (p which, 
as an argument of the wave function, is also an essential quantum degree of freedom. 
The analytic continuation to the Lorentzian wave function 

*(t\<fJ) = yM2H + it\p,f) (2.11) 

takes place by introducing the time-independent Euclidean action I {(f) of the DeSitter in- 
stanton and the Lorentzian Hamilton- Jacobi function S(t,<p), 

I(tt/2H + it, <p) = h{ V ) - iS(t, <p), (2.12) 
I{ V ) = 2I^/2H^)^- w ^- y (2.13) 

and the basis functions of the Lorentzian linearised field equations, 

v n {t) = [u n (n/2H + it)}*, v*(t) = u n (n/2H + it)}, (2.14) 
F^(d/dt)v n (t) = 0, F^(d/dt)5(t - = Sf J^f n (t') ■ (2 ' 15) 

The latter turn out to be the basis functions of the DeSitter-invariant Bunch-Davies vacuum 



24], |32l • The cosmological wave function thus takes the form of the product 

9{t\<p,f) = ^=e-^)/ 2+ ^n^(t,^|/n) (2.16) 



1 Here, for simplicity, we consider the case of a scalar field with F n (d/dT) = — (d/dr)a 3 (d/dr) + while 
for the generic case [^6| a 3 ii n in the exponential should be replaced by the corresponding Wronskian operator 
linear in time derivative acting on u n . 
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of the background wave function and the infinite set of Gaussian states of harmonic oscillators 
representing this vacuum, 

1 



2 



Mt, <p\fn) = -p== exp --SK(t)ti , (2.17) 



n n (t) = -*« 3 (*)||y- (2-18) 



An answer similar to ( |2.1G| ) holds also for the case of the tunneling wave function - the 
only difference is the opposite sign of the Euclidean action in the tree-level part. In what 
follows we shall mainly consider the no-boundary wavefunction, but the main conclusions 
and algorithms for the decoherence aspects of the system will be equally applicable in the 
tunneling case as wel!0. 

An important property of these vacuum states is that their norm is conserved along any 
semiclassical solution 



(Vn.V'n) = J dfn\Mfn)\ 2 = V^A^)]- 1 / 2 , (2.19) 

A n ((p) = ia 3 (v^i) n — Vn v n) = constant . (2.20) 



Note that A n (<p) is just the (constant) Wronskian corresponding to (|2.15|) . We want to 
emphasise, however, that A n is a nontrivial function of the background variable (p, since it 
is defined on full configuration space and not only along a semiclassical trajectory. It can 
thus not be factored out by a deliberate normalisation of these states to unity - a procedure 



that is sometimes assumed in the literature \W\. The replacement of ip n (t, (p\f n ) with the 
normalised states 

"A n \ 1/4 , „ ,, N /ReOA 1/4 (v n \ 1/2 ( 1 



M*Mfn) = (=£) Mt^\fn)={=) (£J exp(-ifi n / n 2 ), (2.21) 
= 1 (2.22) 



would be inconsistent from the viewpoint of the full Wheeler-DeWitt equation or the path- 
integral representation of its solution, because the multiplication of this solution by a con- 
figuration space-dependent object would violate the Wheeler-DeWitt equation. We shall 

2 In the case of the tunneling wave function, the Euclidean path integral derivation does not work, and 
it can be obtained as an alternative solution of the Wheeler-DeWitt equation linearly independent from 
the no-boundary one. Strictly speaking, the behaviours of both of these wave functions in the Euclidean 
and Lorentzian domains does not reduce to a simple analytic continuation presented above. For a no- 
boundary state one underbarrier branch matches with two complex-conjugated branches in the Lorentzian 
domain, while for the tunneling state the situation is reversed. For our presentation here these subtleties are 
immaterial. We shall briefly get back to this issue in the concluding section when discussing the decoherence 
between the branches. 
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give in the rest of this section the justification of the use of the unnormalised vacuum wave 
functions in Q2.17 ). 



2.1. Justification of wave function prefactors: reduction method 
for functional determinants 

One type of justification comes from the general semiclassical theory of constrained 
dynamics [[27| including the Wheeler-DeWitt equation as a particular case. As shown in 



, there exists a unitary map between the solutions of the quantum Dirac constraints and 
physical wave functions in reduced phase-space quantisation. These wave functions, like 
Q2.8|) , have the form of a semiclassical packet with a prefactor given by the Jacobian of the 
transformation from the initial coordinates of the flow of classical trajectories to their final 
coordinates. The matrix of this Jacobian is the solution of linearised equations of motion, 
just like the basis functions in the prefactors of ( |2.8|) , and these basis functions are not 
normalised to unity in the sense of a Wronskian inner product ( |2.20| ). 

A more useful justification of this fact comes from the Euclidean path-integral represen- 
tation of the no-boundary wave function, which goes as follows. The Hartle-Hawking path 
integral in the one-loop approximation reads 

^ 1 E - loop (T\<pJ)= [ D(j)(x)e- I[9{x)] = e~ I{T ^ f) [DetF]" 1/2 . (2.23) 



one— loop 

It contains the tree-level classical action which when expanded up to quadratic order in 
microscopic variables, 

I{r\ V J) = I{rM + \a\r) i ^f n , (2.24) 

gives rise to exponentials of Gaussian oscillator functions; it also contains the one-loop 
functional determinant of the operator 

d(p(x)d(j)(y) 

subject to Dirichlet conditions at the boundary of the Euclidean spacetime ball of "radius" 
t. These boundary conditions follow from the fact that the integration in the path integral 
requires fixed values of fields at the boundary as prescribed by the arguments of the wave 
function. The definition of the functional determinant can then be given in terms of the 
eigenvalues of the corresponding eigenvalue problem 

FMt') = KMr'), < t' < r, <f> a [r) = 0, \ a > 0, (2.26) 
In Det F = Tr In F = ln V ( 2 - 27 ) 
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Here a is a collective index of arbitrary nature enumerating these eigenvalues. 

Suppose now that we can write down an equation which determines these eigenvalues in 
the form 

E(X) = 0. (2.28) 

Its left-hand side E(X) can be used to represent the sum over the eigenvalues as an integral 
in the complex plane of A over the contour C + surrounding the spectrum of the operator - 
so that the residues in the roots of ( |2.28|) reproduce the contributions of different eigenvalues 



Y In X a = — { dX In A 4- In E(X). (2.29) 
„ 2ni J dX 

c+ 

In what follows we assume that this spectrum is positive and lies on the positive real axis. 



Omitting some technical details p8| we integrate in this integral by parts and rotate the 
contour of integration into the contour C_ so that it begins surrounding the origin and 
running over the upper and lower shores of the negative real axis. In view of the absence of 
other roots in the complex plane of A, the only arising residue of the resulting integral gives 
the final answer in terms of the left-hand side of the eigenvalue equation 

£lnA a = --^ /^ln£(A)=ln£(0). (2.30) 
„ 2m J A 

C- 

In the basis of orthonormal spatial harmonics the operator has a diagonal form, so that 
the total equation on its eigenvalues decomposes into an infinite product of partial equations 

F = diag{F n (d/dr)}, E(X) = Y[E n (X), (2.31) 

n 

and the final answer reads 

oo 

£ln A Q = ^lnE n (0) = £ dim(n) ln£ n (0), (2.32) 

a n n =0 

where the summation over the collective index of these harmonics is rewritten as a sum 
involving the degeneracies of the corresponding quantum numbers. 

Let us now get back to the calculation of a particular functional determinant in the 
no-boundary wave function. The role of eigenfunctions in this case is played by the basis 
functions of the operator modified by the extra (mass) term, satisfying the condition of 
regularity at r = 0, 

(F n (d/dr')-X)u niX (T')=0, 0<r' < r. (2.33) 
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Therefore the left-hand side of the eigenvalue equation takes the form 

E n {X) = u nA {r), (2.34) 

and its value at zero argument coincides with the original regular basis function E n (0) = 
v>n,o{T~) = u n {r) taken at r - the "radial" coordinate of the boundary. The use of the equation 
Q2.32|) then immediately gives 

[DetFf^^expf-VlnA^^i^n^^, (2-35) 



U n (T 



where we have also included the contribution of the inflaton sector of the total diagonal 
operator F = diag{F v> (d/ dr) , F n (d/ dr)} . 

The above derivation is presented in a simplified form, because it disregards special 
measures to provide the vanishing of the arc integrals at the infinity of the complex plane 
in the transition from ( 2.29 ) to (|2.30| ), regularisation of ultraviolet divergences and the 



precautions against spurious roots of (|2.32|) caused by multiplication of E n (X) by functions 



of A. These problems are considered in much detail in the authors' paper P^j , where it is 
shown, in particular, that the normalization of the Euclidean basis functions in fl2.34j ) and 
(|2.35 ) should be chosen in such a way that the coefficient of their power-like asymptotic 



behaviour at r = 0, 

w n (r)~[/ n T^, r^O, (2.36) 

is a trivial field- independent constant. Such a normalisation is very different from the nor- 
malisation to unity with respect to the Wronskian inner product. 

3. Wave function norm, Euclidean effective action, and 
the diagonal of the density matrix 

The nontrivial normalisation of vacuum states results in a nontrivial quantum distribu- 
tion of cosmological parameters in inflaton space, yielding a probability distribution for the 



inflaton. This mechanism was used in [18, 29] to analyse the normalisability of the cosmo- 



logical wave function and to find initial conditions for inflation [[5], f|, [19]]. We briefly repeat 
these results here, for they will play a substantial role also in the decoherence properties of 
the density matrix. 

Consider the quantum distribution function of the inflaton field which is just the diagonal 
element of the reduced density matrix obtained by tracing out the microscopic variables, 

p(t\<p)=p(t\<p,<p)= [df\y(t\tpj)\ 2 . (3.1) 



10 



It involves the infinite product of norms (|2.19|) - (|2.20|) that can be reinterpreted by using 
the reduction technique for functional determinants of the previous section. For this, let us 
introduce two sets of Euclidean basis functions inhabiting the total spherical manifold of the 
DeSitter-instanton (the full sphere with the scale factor (|2.5| ) extended to the full range of 
the latitude angle < t/H < ir), 

F E (d/dr)u ± (r) = 0, u~{t) = u(r), w + (r) = u(n/H - r). (3.2) 

They are related by mirror map with respect to the equatorial section of the instanton 
t = n/2H at which the nucleation into the Lorentzian domain takes place. Generically, 
■U-(t) is regular at r = but singular at the antipodial pole of the sphere, r = tt/H, 
and vice versa for u+(t). In the Lorentzian domain they give rise to negative and positive 
frequency basis functions 



u 



[ix/2H + it) = v*(t), u + (tt/2H + it) = v(t), (3.3) 



the Wronskian norm of which can thus be rewritten as a Wronskian product of two different 
Euclidean basis functions calculated at r = ir/2H, 

AM = za 3 (,> n - v* n v n ) t=Q = a\u+^- - ^^)_ /2H • (3-4) 
Let us now modify these functions by the A-term in analogy with the previous section, 

(F n (d/dr) - X)uf(r) = 0, < r < w/H, (3.5) 

still demanding the regularity of u^ij) at r = and of u^{t) at r = tc/H, and define the 
function 

£«M = ^^-^)„ < 3 ' 6 > 
Equation ( 2 . 34j) with this function will be just the requirement of the functional linear 
dependence of (r) and u~^(t) (in view of the second-order equation of motion the equality 
of the functions and their derivatives at one point implies their equality at any r). This 
linear dependence means that there are solutions of ( |3.5|) that are regular on the whole 
instanton including both its poles at particular values of A satisfying this equation. These 
values represent the spectrum of the operator on the whole DeSitter instanton, and according 
to ( |2.32| ) the product of norms arising in ( |3.1| ) reads, up to a numerical constant, in terms 
of this spectrum, 



n(V>n,V>n) = C^exp ^-l^lnA^ = C^exp (-±Tr \nF 



(3.7) 
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It can thus be expressed as the one-loop Euclidean effective action of the theory calculated 
on the DeSitter-instanton of the radius 1/H(cp) |L8|, [19], p9fl , 



A-ioo P (^) = ^TrlnF. (3.8) 

(The total one-loop action contains the contribution of the inflaton mode missing in this 
equation; this explains the origin of an additional A v factor in ( |3.7| ).) 
Thus, the distribution function takes the form 



p(t\<p, <p) = C- — — exp(-I(<p) - ri_i oop (^)). (3.9) 

The derivation of this result was presented here for the case of the no-boundary wave function. 
It is valid, however, also in the tunneling case of the model with microscopic oscillator modes 
The basic difference in this case is the opposite sign of the tree-level Euclidean action 
I(<p) in the exponential. 

The advantage of this representation is that there exist powerful methods of calculation 
and covariant renormalisation of r i-i op( < ^)- These methods, in particular, lead to the 
asymptotic scaling behaviour 

r 1 _ loop (^) = Zln^M, H{<p) ^oo, (3.10) 

with an easily calculable coefficient Z 0, 0. For big positive Z the one-loop contribution 
guarantees the normalisability of the cosmological wave function at large (p. In the model 
with large negative nonminimal coupling of the inflaton it leads to a sharp probability peak 
with the energy-scale parameters belonging to the GUT domain rather than the Planck 
domain 0, |5| and serves as a "source" for inflation. 

4. The density matrix 

We now proceed to discuss the off-diagonal elements of the reduced density matrix that 
are needed for the discussion of decoherence. We shall restrict ourselves to one semiclassical 
branch, i.e. to a wave function of the form (|2.16Q . We shall comment on decoherence between 



branches in section 7. 

The off-diagonal elements of the reduced density matrix read 

p(t\<p, ip') = [ df*(t\<p, /M#', /) , (4.1) 
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where \I/ is given by (|2.17|) . After the Gaussian integration one obtains 



A w A' 



p{t\<p, iff) = C exp [-!/-!/' + i(S - S')} I] [K<{Vn + K)} ^'^ (4-2) 



-1/2 



It is useful to extract from this density matrix the factor y p(t\<p, <p)p(t\(p', <p'), cf. (3.10), the 
remaining factor describing the dynamical entanglement during the Lorentzian evolution. 
By using the relation A n = 2|t> n | 2 Refi n and the expression for the one-loop effective action 
obtained above, we get 

P(t\<p, </) = C /"^^ exp t-\r - \r' + i(S - S')) £>(#, yO, (4.3) 

/V v {t)v'Jt) \ Z Z J 



where 

r = /(^) + r 1 _ loop (^) (4.4) 

is the full Euclidean effective action including the classical part, and 

is the conventional decoherence factor. The form (4.3) of the density matrix corresponds 
to an "ideal measurement" by the environmental variables, since the diagonal part (the 
probabilities) remains unchanged; only the coherence is affected. 

It should be emphasised here that, although the decoherence factor coincides with the 
one known in the literature jTJ], TE, [T^|, the total density matrix ( |4.3| ) is different, because 



it contains the effective- action factors related to both arguments of p(t\tp,tp'). This is the 
effect of a correct normalisation of the vacuum Gaussian states for microscopic variables. 
The factors are entirely defined on the DeSitter-instanton and, therefore, independent of the 
Lorentzian time. Therefore they play an essential role only at the onset of inflation, provided 
the decoherence effects rapidly grow during the Lorentzian evolution. As we shall see, this 
will indeed be the case for all fields except for the case of massless conformally invariant 
fields. 

The amplitude of the decoherence factor can be rewritten in the form 

, Al 1 4Refi„Refi'* . , 

| g (W)l=exp-X> ^f^f , («) 



and the convergence of this series is far from being guaranteed [fl4| , [15| , [16). Moreover, the 
divergences might be not renormalisable by local counterterms in the bare quantised action. 
We shall now analyse this question of coherence for all types of bosonic fields; the case of 



fermions will be dealt with in a separate paper El] . 
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5. Inconsistencies in the dimensionally regularised den- 
sity matrix 

In this section we shall show that dimensional regularisation leads to density matrices 
that violate crucial properties and that are therefore inconsistent. We shall start with the 
case of a massive minimally coupled scalar field. There the equation for the Lorentzian basis 
functions reads 

= 0. (5.1) 

In what follows we shall use dimensional regularisation to regulate the divergencies. In this 
regularisation the equation retains its form except for the power of a in the first term, which 
instead of 3 becomes d — 1 where d — > 4 is the regularisation parameter of the spacetime 
dimension (and n is being changed to n + (d — 4)/2). The corresponding solutions for 
the DeSitter background, with the Hubble constant H, that are regular on the Euclidean 




hemisphere of the gravitational instanton (at r = n/2H + it = 0) read in this dimension |28 



v ni t) = (cosh i/ t )-^p : ^^__( i si „ h my (5.2) 



The behaviour of the series ( f4.6|) depends on the behaviour of fl n with respect to n. Since 
we shall be strongly interested in the limit of large masses, it is worth presenting the above 
function in the form of an 1/m asymptotic expansion uniform in n p8| , [31]. It reads 

- 1/2+iAl J |r(n + iA + l/2)| VX 

x exp ( — In Aarctan— ^ fc ^ ,' — -, A — > oo, (5.3) 

V2 1 + v a) ^ A fc ' ' y J 

where 

a =yP = l + a,v= (p 2 _ z 2y/2 > ( 5 - 4 ) 

and the coefficients Tj.(v,a) beginning with T (v,a) = 1 are given in (they are not 
needed below). They are uniformly bounded for all values of A and /i. In view of this 
boundedness, the exponential part of this function represents for large n the conventional 
high-frequency adiabatic expansion and contributes the dominant part to Q n , 



fL = a 



Vn 2 + m?a 2 + i sinh Ht 1 



1 ma 



2„2 



2 n 2 + m 2 a? 



+ 0(l/m). (5.5) 
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The corresponding leading contribution to the amplitude of the decoherence factor, 



1 ^ 4a Wn 2 + m 2 a 2 v^T 



2„/2 



In \D(t\<p, - T E " ln 7 , ~ , (5-6) 



4 



n=o (a 2 \/n 2 + m 2 a 2 + a! 2 \Jn 2 + m 2 a' 2 



obviously contains cubic and linear divergences which cannot be represented as additive 
functions of a and a'. This means that no one-argument counterterm to r and r' in ( |4.3|) 
can cancel these divergences of the amplitude - an observation made a number of years ago 
in 



Let us forget for the moment this difficulty and just use the regularisation which iden- 
tically discards all power divergences. Then, provided that the logarithmic divergences are 
vanishing, one can at least have a finite quantity subject to physical analysis. Another prob- 
lem can, however, occur in this case - this regularisation can perform an "oversubtraction" 
in the sense that the resulting regularised density matrix will be inconsistent - it will vio- 



late the necessary properties satisfied for a convergent series (4J5). The calculations briefly 
presented below really show that the dimensionally-regularised density matrix would have 
off-diagonal terms infinitely growing for ip — <p' — > oo and thus being inconsistent; for ex- 
ample, Tr p 2 would diverge. We emphasise that this problem has not been discussed before, 
since reduced density matrices are usually not considered in quantum field theory. 

Dimensional regularisation of the divergent sum ( |5.6p implies that it should be reformu- 
lated in d-dimensional spacetime, which boils down to two main modifications - the degener- 
acy of the n-th eigenvalue of the spatial Laplacian n 2 , and the eigenvalues A n = —n 2 /a 2 itself 
should be replaced by their d-dimensional versions dim(n, d) and X n (d). From the theory of 
irreducible representations of the group 0(d + 1) it is well known that 

(2n + d-4)T(n + d-3) 



dim(n, d) 



T{n)T(d-l) 



„ t {2-e/n)(l-e/n)T{n-e)n* ,_ £ 



n rW \ r( \ n ' e = 4 " rf ' ( 5 - 7 ) 

T(3 - e) r(n) 



while the d-dimensional A n (d) can be obtained from A n by shifting the quantum number n, 
n — > n — (4 — d)/2. The latter results in the following modification of the frequency, 



^n(d) = a 



2 

\ 



4-<f 2 



n — +m 2 a 2 + 0(m). (5. 



Note that the main effect of regularisation in ( |5.7| ) consists in the multiplication of the 
summation measure n 2 by the factor n~ e which for large positive e = 4 — d can make the 
divergent series convergent. The analytic continuation of the result to e = allows one to 
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disentangle the divergent part of the whole answer as a pole in dimensionality and get a 
finite contribution. 

The further calculations are based on the use of the summation Euler-Maclaurin formula, 



E/( 



n 



n=l 



dyf{y)-^f{0) + i 



dy 



f(w) - f(-iy) 



(5.9) 



e 2ny _ I 

which allows one to reduce the infinite sum to the combination of integrals along real and 
imaginary axes in the complex plane of n = y. Because of the exponential damping factor in 
the integral over imaginary axes, the regularisation can be removed in its integrand from the 
very beginning. This would mean that its contribution is vanishing (because the integrand is 
an even function of its argument), except for the possible residues at the singular points of the 
integrand given by the equation n 2 + (ma) 2 = —y 2 + (ma) 2 = 0. These residues will, however, 
contribute the exponentially small terms 0(e~~ ma ) negligible in the approximation of large 
masses ma ^> 1. Thus, the decoherence factor reduces in the main to the contribution of 
the integral over the real axes, which after the change of the integration variable y = ma/x, 
takes the following form, 

\n\D\= dxx e - 4 f(x,e) + 0{e~ ma ). (5.10) 
Jo 

Its integrand can be obtained from (|5.6| ) by the substitution of n = ma/x and the dimensional 
modifications of the above type, and for nonvanishing e it looks rather complicated. However, 
its complexity due to the nontrivial dimensionality can be represented in the following form, 
2(ma)~' 



f(x,0) 



r(3 
i 



-f(x, 0) [l + 0(e/ma) 



4 



ma) In • 



48 2 VTT^ 2 VI + S 2 x 2 



2 i 



(5.11) 
(5.12) 

by introducing the new parameter 5. 



S 



a 
a 



'VTTx 2 + 8 Vi + 5 2 x 2 ] 

with the function f(x, 0) trivially obtained from 
Due to this, the resulting expression does not look symmetric in a and a' , but it will allow 
us easily to analyse the behaviour of far off-diagonal elements corresponding to 5 <C 1. 

The further transformations of ( |5.10p consist in the consecutive integrations by parts in 
the convergence domain of the integral and the analytic continuation to e = 0: 

f (0,0) 



dxx e f(x,e) 



d 



1 r°° 

-I dx In x / 

6 jo 

f"(0, e) 



(iv), 



x, 0) 



de ( e -3)(e-2)(e-r 



e=0 



In view of ( |5.11| )-( |5.12| ) one has 



no,o)=o, 



S A0,e) 



{ma) 3 0{l/m), 



(5.13) 



(5.14) 



e=0 
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which means the absence of logarithmic divergences, and the final result reads 

roo d 4 f(x 0) 7T 

ln|D|= / dxlnx , ' = — (ma) 3 + (m 2 ) , a > a'. (5.15) 
Jo dx A 24 

This implies the inconsistency of the regularised density matrix in view of its infinitely 

growing off-diagonal elements. 

The behaviour of the regularised density-matrix elements in close vicinity of its diagonal 

is also unsatisfactory. For small \tp — (p'\ the dimensionally regularised decoherence factor in 

the limit of large mass gives the result 

In \D(t\<p, <p')\ ~ ^m 3 d(a - a') 2 , (5.16) 
o4 

also leading to the inconsistency of the density matrix. 

The physical inconsistency of the dimensional regularisation is even more transparent 
in applications to massless conformal invariant fields. On Friedmann backgrounds, such 
fields decouple from the gravitational variables, so one should expect the complete absence 
of decoherence effects, D(t\ip,(p') = 1. This is, however, not the case. Simple calculations 
show that for a massless conformally coupled scalar field the basis function v*(t) and the 
frequency function have the form: 

n 

^^b^hi ' (5 - 17) 



tt n = na 2 + ia 2 VH 2 a 2 - 1. (5.18) 
The corresponding decoherence factor is then given by the divergent sum, 

ln|D| = --£nMn A + - , (5.19) 



4 ^ V n 2 

V 2aa' J ' V 2aa> J 

Its dimensional regularisation gives the exact result 

ln|D(%,yy)| = ^(£) 3/2 , (5.21) 

again demonstrating the wrong behaviour of the off-diagonal elements of the density matrix 
in view of the obvious growth of the exponential (B/A) 3 ^ 2 — > oo, \ip — (p'\ — > oo. For a 
minimally coupled scalar field the situation is analogous to the conformal case. 

Thus, in the massless limit the regularised decoherence matrix turns out to be as in- 
consistent as for massive fields. Moreover, for conformally invariant fields it leads to an 
anti-intuitive conclusion of growing quantum correlations between the environment and grav- 
itational background in spite of their actual decoupling. 
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6. Conformal parametrisation of bosonic fields 



In this section we discuss a possible solution to the problems encountered with di- 
mensional regularisation. It was remarked in the Introduction that a redefinition of the 
environmental fields can change the reduced density matrix. Moreover, as was remarked in 



17|| ) such a redefinition can be motivated by physical considerations. In fact, it was already 



shown for a specific model in [O] that a redefinition can lead both to a finite decoherence 
factor and to a consistent density matrix. Here we shall demonstrate that such a viable 
redefinition is not arbitrary and can lead to a unique result. 

Let us consider the following redefinition of a bosonic scalar field, 



f{t) - f(t) = a»(t) f(t) 



(6.1) 



with an arbitrary numerical parameter /i. Under this redefinition the kinetic term of the 
Lagrangian changes as 



L(f) 



'df 
dt 



+ 



7 3-2/i 



df 

dt 



7x 2 



+ 



and the Wronskian operator and basis functions get replaced by 



W(d/dt) = a: 



d_ 

dt 



W{d/dt) = a 



3-2/1 



d 



v n (t) v n (t) = a^(t)v n (t). 
The new frequency function then reads 

tin = -ia^Unv* n 
dt 



•2) 



(6.3) 
(6.4) 



,2-2/. 



1 



m 2 a 2 



Vn 2 + m 2 a? + i sinh Ht [ 1 — A* H - - 

\ 2 n 1 + m 2 a 2 



+ 0(l/m) 



(6.5) 



For a particular choice of the parameter [i = 1 two important things happen with this 
function - the leading n behaviour of Cl n becomes independent of the macroscopic variable, 
and its imaginary part decreases with growing n as l/n 2 , 



Reb n 
Im b r , 



Vn 2 + m 2 a? + b n 



0{l/n) 



sinh Ht 



m 2 a 2 



+ 0(l/m) = 0{l/n 2 ). 



2 n 2 + m?a 2 

Substituting this behaviour into ( f4.5|) , one finds 



ln£>(%V)^ 7 E 



.Im&l 



. Im b Tl 



n 



— i- 



n=0 



n 



+ 0(l/n 4 ) 



(6.6) 
(6.7) 

(6.8) 
(6.9) 
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This expression has at most logarithmic divergences which are imaginary and, thus, affect 
only the phase of the density matrix. Moreover, since b n and b' n depend on a and a', respec- 
tively, these divergences decompose into an additive sum of one-argument functions and, 
therefore, can be cancelled by counterterms to the classical action S (and 5") in ( |4.2|) . The 
real part is simply convergent and generates the finite decoherence amplitude. This result is 
formally similar to the result for the decoherence factor in QED [15| . 

Let us show now that no inconsistencies arise with this redefined density matrix. To 
begin with, let us consider the far off-diagonal behaviour of the decoherence-factor amplitude 
for the arbitrary parameter of reparametrisation \i in the large-mass limit. In this limit 
£l n ~ a 2 ~ 2 ^\/n 2 + m 2 a 2 , and the decoherence factor can approximately be represented by 
the integral (replacing the summation, n = may) 

4<5 2 -^ v /TTFv / 2/ 7 T7 2 ~ 



\n\D(t\<p,<p') 



[ma) 







dyy In 



VT+y 2 + 5 2 ~ 2 ^y 2 + 5 2 ) 



2 ■ 



(6.10) 



where 5 = a' /a. This integral is convergent only for fi — 1, and for a ^> a' it yields the 
answer which shows the suppression of far off-diagonal terms, 

irr- 

3, 



\D(t\tp,tp')\ ~ exp 



[ma) 



24 



0(m' 



(6.11) 



For other values of fi this integral should be regulated to give the final answer. Dimensional 
regularisation in the main changes the power of y in the integration measure from y 2 to y 2 ~ £ . 
Then the integration gives two different results for different ranges of /i 



\£>(t\<p,<p') 

\D(t\<p,<p') 



exp 



exp 



7T 

24 



(ma) 



71 



24 



(ma) 



, fi > 1, 
li < 1 . 



(6.12) 
(6.13) 



In the latter case this obviously coincides with the behaviour of the inconsistent dimensionally 
regularised density matrix (|5.15|) that was naively defined for /i = 0. Thus, although a 
sensible result can be obtained for ji > 1, the value \i — 1 is distinguished in yielding 
directly a finite result. In spite of this, the results for /i = 1 and /i > 1 are qualitatively the 
same. 

In the vicinity of the diagonal the modified density matrix also acquires a good quasi- 
Gaussian behaviour. Expansion of the expression given above in powers of small \a — a'\ 
gives for large m the dominant result 

, . ~ , . , XI m 3 7ra(a — a') 2 

In \D(t\<p, <p')\ = K - (6.14) 



64 
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For massless conformally-invariant fields, the modified density matrix justifies the expec- 
tations induced by the decoupling of the environment from the gravitational background, cf. 
[ DJ . Its decoherence factor turns out to be trivial. This can be seen by noting that, say, 



in the case of massless conformally-invariant scalar field, the basis and frequency functions 
have the form: 

1 + zsinhff^ 2 

l-.sinhgt j ' (6 ' 15) 
Q n = —ia-r hi^nW = n ■ (6.16) 

(JjL 

Hence, D(t\ip,ip') = 1. The same holds also for the electromagnetic field. 

We conclude this section by briefly considering two massless conformally non-invariant 
fields - a minimally coupled scalar field and gravitons. What they have in common are the 
basis- and frequency functions in their respective conformal parametrisations: 

rt 

(l + ismhHt\ 2 (n - z sinh if A 
<® = \ l-i^Ht ) { n+l J' (6 - 17) 

They differ only by the range of quantum number n (2 < n for inhomogeneous scalar modes 
and 3 < n for gravitons) and by the degeneracies of the n-th eigenvalue of the Laplacian, 

dim(n) scal = n 2 , (6.19) 
dim(n) grav = 2{n 2 - 4). (6.20) 

The decoherence factor generated by both of these quantum fields thus has the form of the 
following, obviously convergent, series: 

2 



\n.\b{t\tp,ip')\ = — ^dim(n) In 

4 n=2 



( sinh 2 Ht - sinh 2 H't 

1 + 



4(n 2 + sinh 2 Ht)(n 2 + sinh 2 H't) 

(cosh 2 Ht sinh Ht (n 2 + sinh 2 H't) - cosh 2 H't sinh H't (n 2 + sinh 2 Ht) 
An 2 (n 2 - l) 2 (n 2 + sinh 2 Ht)(n 2 + sinh 2 H't) 



.(6.21) 



The calculation of the corresponding far off-diagonal elements is more complicated here, 
but qualitatively it reproduces the behaviour (6.11) for a massive scalar field with the mass 
parameter replaced by the Hubble constant of the gravitational background, 

\D(t\<p,<p')\ ~e- c{Ha) \ a>a', C > 0. (6.22) 
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To study the behaviour of matrix elements in the vicinity of diagonal of density matrix 
it is enough to expand the expression (|6.21 ) up to the second order in the quantity (H — H') 
and calculate the arising series by using []32] . The exact result even in this order of (H — H') 
is too long to be fully presented here. Thus we restrict ourselves to the limit of late time t 
for this expression, 

\D(t\ V , V ')\~exp^(H-H>)H 2 e 4Ht y Ht » 1, (6.23) 
which clearly shows a rapid growth of decoherence during the inflationary evolution. 



7. Conclusions 

Let us recapitulate briefly the main results of the present paper. We have calculated the 
reduced density matrix for the inflaton field in a model of chaotic inflation by tracing out 
inhomogeneous degrees of freedom corresponding to various bosonic fields. In the original 
parametrisation of these fields, ultraviolet divergences arise in the decoherence factor of the 
density matrix. The dimensional regularisation of these divergences was shown to violate the 
consistency of a reduced density matrix as a bounded operator, which is apparently related 
to the fact that these divergences do not have the structure of one-field counterterms usually 
used for ultraviolet renormalisation. A physically motivated conformal redefinition of the 
environmental bosonic fields leads, however, to well-defined finite results, which show that 
the Universe acquires classical properties near the onset of inflation. 

The seemingly mysterious finiteness property for the decoherence factor in the confor- 
mal parametrisation has a natural explanation. The off-diagonal elements of the density 
matrix, in the language of Feynman diagrammes, are roughly proportional to the series of 
gravitational-matter vertices with two quantum matter legs and the growing number of the 
gravitational (background fields) ones. Indeed, when expanded in Taylor series in a — a' (and 
if — if'), the one- loop density matrix can be represented as a quantum loop with the insertion 
of external legs corresponding to the differentiations with respect to background fields, these 
insertions forming the vertices in question. Such vertices generally contain two derivatives 
which make the loop graph strongly divergent. The conformal parametrisation of quantum 
matter fields on the Friedmann background effectively removes the gravitational variables 
from the kinetic term of the matter Lagrangian, thus actually removing the derivatives from 
the corresponding vertex. This makes the relevant one-loop Feynman diagrammes finite. 

The obtained results show that the decoherence properties of the system strongly depend 
on the parametrisation of the quantum fields that are traced out. Although, intuitively, 



21 



this sounds unsatisfactory - the observable correlations depend on the choice of variables 
about which we all the same do not have any information - this situation is very typical in 
quantum theory. In fact, the density matrix is analogous to the S-matrix taken off shell. 
As is well known, quantities such as the off-shell S-matrix, the off-shell effective action, 
strongly depend on the parametrisation of quantum fields, the integration over which is 
performed when calculating these quantities. A similar situation arises here: In contrast to 
the diagonal elements given by the on-shell effective action, the off-diagonal elements are 
ambiguous. For the bosonic fields that we considered here, this situation seems to have 
a relatively satisfactory resolution - there exists a distinguished parametrisation in which 
divergences are absent. This parametrisation is distinguished for several reasons. First, it 
reflects the conformal weight properties of the field in question even in the case when this 
field is not conformally invariant. Secondly, this parametrisation is attained by an overall 
multiplication of the field by some power of the scale factor, which does not spoil the locality 
of the field. As we shall show in a forthcoming paper []21j ], the situation with fermions is 
less pleasant. Conformal parametrisation leaves the answers divergent, while the additional 
reparametrisation that can make them ultraviolet-finite is nonlocal - it has the form of 
Bogoliubov transformations rotating the different harmonics of the field differently, thus 
destroying the local nature of the original field. 

The results obtained have strong implications for the quantum-to-classical transition 
within the theory of the quantum origin of the inflationary Universe. They can be especially 
important in the quantum cosmology of the chaotic inflationary model with big negative 
nonminimal coupling £ of the inflaton field to curvature 0, |l^, || This model was 

shown to generate a probability peak in the distribution function of the inflaton field p - the 
diagonal element of the density matrix - at the GUT energy scale for the inflationary Hubble 
constant H(tp) ~ m-p ~ l(T 5 mp (A is a coupling constant in the quartic term of the 

inflaton potential, and the value of the ratio Va/|£| ~ 10~ 5 ~ AT/T follows in this model 
from the COBE normalisation of the magnitude of the microwave background anisotropy). 
This probability peak has at the onset of inflation a very narrow quantum width for the 
Hubble constant, AH/H ~ >/A/|£|. Now we can show that the width of the density matrix 
in the off-diagonal direction at the initial moment of time t = has the same magnitude. 
This follows from ( |4.3|) . 

First, apart from the decoherence factor D(0\(p, if') the only strongly suppressing factor 
in this equation is, cf. ( |3.10| ), 

ex P (-r 1 _ loop /2 - r;_ loop /2) ~ {hh')- z '\ (7.1) 

with a very big parameter of anomalous scaling in this model, Z ~ |£| 2 /A. It gives a power- 
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like suppression of the far off-diagonal terms, and the quantum width in the off-diagonal 
direction, in the vicinity of the diagonal, 

(HHT Z/2 ~ H- z e -(^')VAH^ (7 2) 

AH ~ Hy\/\£\, coincides with the width of the wave packet. 

As far as the decoherence factor D is concerned, at the start of inflation it is basically 
inefficient, D(0\(p,(p') = 1, for a very simple reason. Large |£| generate in this model by the 
Higgs mechanism very big values of masses for all quantum fields directly coupled to inflaton. 
The order of magnitude of their mass parameters is such that the ratio m 2 /H 2 ~ |£| is very 
big and (p- independent. This means that for big masses at t = the dominant argument 
of the decoherence factor mo(0) = m/H is ^-independent and D(Q\tp,(p') = 1. Thus, the 
Universe is essentially quantum at the "start of inflation" . But due to decoherence effects it 
rapidly becomes classical: the quantum width of the distribution function - the probability 
- stays basically the same, while the off-diagonal width rapidly decreases. From (|6.14|) it 
immediately follows that this width for late time is AH ~ exp(— 3i/t/2)/(t|£| 3 / 4 ). 

In our paper we have restricted ourselves to decoherence within one semiclassical branch 
of the wave function. There is, however, also the possibility that different branches might 
interfere. For the examples discussed in our paper, this would be the case for the no- 
boundary state that is a superposition of two complex conjugate semiclassical branches. 
Decoherence is, however, also effective in suppressing interferences between such branches as 
long as decoherence within one branch holds At the technical level, this follows from the 
expression (|4.5|) which for the interbranch case has in the denominator Q n + Q' n rather than 
Q n + Q'*. This means that the imaginary parts of the frequency functions add up instead 
of partially cancelling one another. Therefore, the amplitude of the interbranch decoherence 
factor is smaller than that of one branch. In quantum mechanics, a nice analogy is provided 
by the case of chiral molecules where the superposition between the left-handed and the 
right-handed version of the corresponding molecule is being suppressed after the localisation 
has been taken place by continuous measurement with, for example, light in each version. 

We have here not included a discussion of the influence of fermionic fields, since this 
will need a different formalism and reveal novel aspects. This case will be discussed in a 



separate paper |fnj that, together with the present paper, should present a complete picture 
of decoherence in one-loop quantum cosmology at the onset of inflation. 
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